function y = lorentzian(p,x)
    y = p(1)*p(2)./((x - p(3)).^2 + p(2).^2) + p(4);
end